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ABSTRACT 

Based on the characteristics of the magnetorotational instability (MRI) and the MRI-driven turbulence, we 
construct a steady model for a geometrically thin disk using "non-standard" a-prescription. The efficiency 
of the angular momentum transport depends on the magnetic Prandtl number, Pm = v/r), where v and r\ are 
the microscopic viscous and magnetic diffusivities. In our disk model, Shakura-Sunyaev's a-parameter has 
a power-law dependence on the magnetic Prandtl number, that is a oc Pm s where 8 is the constant power- 
law index. Adopting Spitzer's microscopic diffusivities, the magnetic Prandtl number becomes a decreasing 
function of the disk radius when 8 > 0. The transport efficiency of the angular momentum and the viscous 
heating rate are thus smaller in the outer part of the disk, while these are impacted by the size of index 8. We 
find that the disk becomes more unstable to the gravitational instability for a larger value of index S. The most 
remarkable feature of our disk model is that the thermal and secular instabilities can grow in its middle part 
even if the radiation pressure is negligibly small in the condition 8 > 2/3. In the realistic disk systems, it would 
be difficult to maintain the steady mass accretion state unless the /^-dependence of MRI-driven turbulence is 
relatively weak. 

Subject headings: accretion, accretion disks — instabilities — magnetic fields — MHD 



1. INTRODUCTION 

The physical mechanism to transport angular momentum in 
accretion disks is an important unsettled issue in astrophysics. 
It has been extensively studied for decades using theoretical 
and numerical procedures. The main purpose of these studies 
is to find the efficient transport process of the angular momen- 
tum to account for the powerful mass accretion deeply asso- 
ciated with the release of th e gravitational en ergy observed in 
a strophysical disk sys tems (Ki ng et alJl2007l) . 

Shakura & Sunvaevl (1 19731) propose a pioneering disk 
model for describing the steady mass accretion, in which the 
efficiency of the angular momentum transport is parametrized 
and represented by a dimensionless viscosity parameter a = 
—trcj, I p, where t r< j, is the r0-component of the stress tensor, 
and p is the total pressure. The phenomenological parame- 
ter a introduced in the disk model should be determined ex- 
perimentally. This disk model is widely accepted today as 
the stan dard, and referred to as "Standard Shakura-Sunyaev 
Model" dKato et al.l [l998. hereafter, we simply call it the 'a- 
model'). 

It is pointed out by King et al. (2007) from the observa- 
tional point of view that the viscosity parameter a should 
be at the level a ~ 0.1 in order to explain the powerful out- 
bursts from disk systems, such as dwarf novae, X-ray tran- 
sients, and a fully ionized part of the disk in Active Galac- 
tic Nucleus (AGNs). Since the required level for the viscos- 
ity parameter a is much larger than that resulting from the 
microscopic molecular viscosity, the turbulent viscosity be- 
comes the most promising candidate for governing the angu- 
lar momentum transport in astrophysica l disk systems (see 
also lHartmann et all 1 19981 ISmakl Il999t iDubus et ail 120011; 
IStarling et al.ll2004HLodato & Clarkdl2004l) . 

The Keplerian disks are known to be linearly stable to the 
shear instability (Rayleigh instability) which could power the 
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hydrodynamic turbulence. Even when the non-linear per- 
turbations are imposed, the Keplerian disks can never reach 
the highly turbulen t state in the Rayleigh stable condition 
(Balb uset al.lll996l) . Numerical studies in addition show that 
the other hydrodynamic instabilities, such as convective and 
Papaloizou-Pringle instabilities, do not play a central role 
in inducing th e efficient angular momentum transport in the 
disk sy stems dPapaloizou & Pringlel Il984t iBlaes & Hawleyl 
[19881 IStone & Balbusl Il996t Balbus & Hawlevl 11998 : 
Ijqhnson & Ga mmie 200^). 

Bal bus & Hawlevl d 1 99 lb focus, for the first time, on 
the mag neto-rotational i nsta bility (MRI) wh i ch is originally 
found by Velikovl (ll959l) and [Chandrasekhai] dl960|) as a can- 
didate for operating the magnetohydrodynamic (MHD) turbu- 
lence in the disk systems. It has essentially a local nature and 
destabilizes weakly magnetized differentially rotating sys- 
tems with a negative shear rate q = dlnfi/dlnr, where f2 is 
the angular velocity and r is the cylindrical radius. The non- 
linear MHD turbulence driven by the MRI has been broadly 
investigated as the leading mechanism for the turbulent an- 
gular momentum transport in disk systems using state-of- 
the-art numerical tech niques and huge computer facilities for 
the MHD simulations dHawlev & Balbus E992t lHawleyl 19951 
ISano et aDl2004t IFromang et al.ll2007l) . The physical proper- 
ties of the MRI-driven turbulence, however, have not yet been 
confirmed completely even today. 

A recent remarkable finding in the MRI studies is 
the dependence of the transport efficiency of the MRI- 
driven turbulence on th e ma g netic Prandtl n umber 
Pm dLesur & Longarettil 120071: IFromang et all 120071: 
I Simon & Hawleyl 1200 9). where Pm = v/r), v is the vis- 
cosity and r] is the magnetic diffusivity. The local shearing 
box simulations explicitly taking into account the micro- 
scopic kinematic viscosity and resistivity suggest that the 
total turbulent stress (sum of Maxwell and Reynolds stresses) 
depends on the magnetic Prandtl number. Especially in the 
system with the non-zero net magnetic flux, the turbulent 
stress and the viscosity parameter a resulting from the MRI- 
driven turbulence would increase with the magnetic Prandtl 
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number, that is a oc Pm s , where <5 is the power-law index 
with positive value (Froman g et al.l l20"07t ISimon & Hawleyl 
l2009h . 

It should be stressed that the viscous effect on the non-linear 
MRI-turbulence is seemi ngly contradictory to that on the lin- 
ear gr owth of the MRI (iMasada & Sanol2008tlPessah & Chaiil 
2008). Surprisingly, the viscosity can boost the transport ef- 
ficiency of the MRI-driven turbulence although it suppresses 
the growth of the MRI in its linear evolutionary regime. The 
physical mechanism for the viscous boosting of the MRI- 
driven turbulence is a central issue in the current MRI studies 
and remains to be solved to draw a complete physical pic- 
ture of the angular momen tum transport in the accretion disks 
dLongaretfi & Lesudl20lol) . 

The impact of the microscopic diffusivities on the MRI- 
driven turbulence should be essential not only in understand- 
ing the nature of the MRI itself, but also in applying it to 
the realistic astrophysical disk systems. Assuming the stan- 
dard a-disk, the size of the magnetic Prandt l number drasti- 
cally c hanges depending on the disk radius. Bal bus & Henril 
(2008) suggest that the magnetic Prandtl number Pm is typi- 
cally larger in the inner disk when we adopt Spitzer's value s 
for evaluating the microscopic diffusivities (Spitzer 1962). 
Hence, the transport efficiency of the angular momentum 
should be affected by the magnetic Prandtl number and be- 
come a function of the disk radius when the viscosity pa- 
rameter follows the relation a oc Pm s w hich is implied from 
the recent local shearing box s imulations dLesur & Lon garetti 
l2007t ISimon & Hawlevll2009h . 

In this paper, we construct a steady model for geometri- 
cally thin disk taking account of the "non-standard" descrip- 
tion of the viscosity parameter which depends on the magnetic 
Prandtl number in the form a oc Pm 5 . This should be consis- 
tent treatment of the turbulent viscosity with the results from 
recen t local shearing box si mulations (Lesur & Longaretti 
l2007t ISimon & H awlev 2009). A central difference between 
the work done by iBalbus & Henril (120081) and ours is to in- 
clude the structural change of the disk depending on the size 
of the magnetic Prandtl number Pm in our model consistently. 

This paper is organized as follows: The steady disk model 
with the non-standard description of the viscosity parameter 
a oc Pm b is constructed in § 2. We then describe the basic 
properties of our disk model. In § 3, we investigate the stabil- 
ity of our steady disk model to the gravitational, thermal, and 
secular instabilities. The dependence of the stability criteria 
on the magnetic Prandtl number (or index <5) is our interest in 
this section. Finally, we discuss the application of our model 
to the realistic astrophysical disk system and summarize our 
findings in § 4. 

2. GEOMETRICALLY THIN DISK WITH NON-STANDARD 
q-PRESCRIPTION 

We construct a steady disk model consistently com- 
bining the properties of the M RI-driven turbulence re- 
cently found by numeric al studies (Lesur & Lon garettill2007l: 
Simon ~& Hawlevl 12009 ^ . The fundamental assumption 
adopted in this work is that the viscosity parameter a de- 
pends on the disk radius and reflects the transport properties 
of the MRI-driven turbulence. Also we assumed that the disk 
heating is local in which the variation of the thermal energy 
responds to the turbulent en ergy dissipation instantaneously 
(Balbus & Papaloizoul 1 1 9991) . This assumption is consistent 
with the recent numerical work of the MRI-driven turbulenc e 
dSano & Inutsu kal200lHSimon et al]|2009HGuan etal]|2009l) . 



According to the recent numerical studies, the efficiency of 
the angular momentum transport sustained by the MRI-driven 
turbulence is regulated by the magnetic Prandtl number which 
characterizes the system. Using the viscosity parameter a, 
the transport efficiency of the MRI-driven turbulence can be 
represented by, as is described in § 1, 

a = aoPm s , ( 1 ) 

dLesur & Longarettil2007l:ISimon & Hawlevll2009l) . where a 
is the normalization parameter which controls the size of a to 
be smaller than unity. Note that this relation reduces to the 
classical standard a-model by taking 6 = 0. 

For providing the magnetic Prandtl number, we adopt the 
Spitzers' values for the microscopic kinematic viscosity v and 
the electron resistivity r\ by assuming the fully ionized gas. 
This assumption is reasonable for the disk of the X-ray bina- 
ries and for the inner part of the disk in AGNs. Then v and r\ 
are expressed as 



and 



i/=1.6x 10 -I V TH\nA HH ) 
?7 = 5.55 x 10 n HlnA e// , 



(2) 



(3) 



where p is the disk mass density, and T is the disk temper- 
ature (Spitzer 1962). The Coulomb logarithms for proton- 
proton and electron-proton scatterings are represented above 
by lnA//H and lnA^, respectively. The magnetic Prandtl 
number is then described as 



Pm = 



4.2 x 10 4 K 



10 14 crrT 
nl 



(4) 



where n is the number density and I = lnA e //lnA#H is the 
product of the Co ulomb logarithms and fixed to be I = 40 in 
the following (see lBalbus & H enri 2008|). 

The governing equations are the conventional used in the 
classical standard a-model: 



Mass conservation equation 

M = -2?rrv r E . 

Force balancing equation (Keplerian rotation) 



n = n K = 



GM 



Angular momentum conservation 




Energy balancing equation 

9 2 32ar 4 

— VfiAr = . 

4 3r 

Hydrostatic balance in the vertical direction 

n 



(5) 



(6) 



(7) 



(8) 



(9) 



where J7k is the Keplerian angular velocity. The physical pa- 
rameters M, M, w, c„ v t , H, R m and r are the mass of the 
central objects, the mass accretion rate, the radial velocity, the 
sound speed, the turbulent viscosity, the vertical scale height 
of the disk, the inner radius of the disk, and the optical depth 
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in the vertical direction. Here S is the surface density defined 
by £ = 2pH. The physical constants G, c, and a are the gravi- 
tational constant, the speed of light, and the Stefan-Boltzmann 
constant respectively. 

In this paper, we naively assume that the viscosity parame- 
ter a is related to the turbulent viscosity by a relation, 



log Pm 



(10) 



(but, see ISano et a"ill2004t IPessah et al.ll2007l and §g]in this 
paper). Note that there is a clear distinction between the turbu- 
lent viscosity v, and the kinematic viscosity v. In our model, 
the turbulent viscosity v t is determined by the MRI-driven tur- 
bulence which controls the accretion dynamics. The transport 
efficiency due to the MRI-driven turbulence is regulated via 
the relation (fTJ by the microscopic kinematic viscosity v. The 
kinematic viscosity affects the global nature indirectly by reg- 
ulating the MRI-driven turbulence. 

The opacity of the disk k, which relates to the optical depth 
rasK = 2t/£, is given by 

= K es + K pT~i . (11) 

where K es ~ 0.4 cm 2 g" 1 is the opacity due to the electron scat- 
tering and Kff is that contributed from the free-free absorption 
with kq ~ 6.4 x 10 22 in cgs unit. 

For the closure of the system, we need the equation of state 
for the composite gas of the photon and baryon, 

4aT 4 



P 
fim p 



k B T+- 



3c 



(12) 



where m p is the proton mass, ks is the Boltzmann constant and 
[i is the mean molecular weight providing /i = 0.5 throughout 
this paper. 

Solving coupled equations (|T), (l4l- (fT2l i with the zero 
stress boundary conditio n at r = Ri n (see, equation [7] or, 
iNovikov & Thornell 1973b . we can obtain the local disk struc- 
tures for three characteristic regimes. In the inner part of the 
disk, the radiation pressure dominates the gas pressure. In 
contrast, the gas pressure becomes dominant in the middle 
and outer regions of the disk in our model. The difference be- 
tween the middle and outer regions is the source responsible 
for the opacity. The electron scattering plays a crucial role in 
determining the opacity in the middle region. The outer struc- 
ture of the disk depends mainly on the opacity contributed 
by the free-free absorptio n (see, also standard a-model in 
Sha kura & Sunvaevlll973l) . 

The steady disk structure obtained here is then character- 
ized by four dimensionless parameters, x = r / r s , m = M/Mq, 
m = Mc 2 /L E and 5 = log(a/a )/logPm, where r s is the 
Schwarzschild radius, M Q is the solar mass, and L E is the 
Eddington luminosity, respectively. The innermost radius R m 
is assumed as R[ n = 3r s in the following. 

Figures [T] and [2] show the radial profiles of the magnetic 
Prandtl number (left top), a (right top), the radial velocity 
in unit of c (left middle), the scale height (right middle), the 
surface density (left bottom), and the disk temperature (right 
bottom) for the models m = 10 (Fig. Q] the model X-ray bi- 
naries ) and m = 10 8 (Fig. [2] the model AGNs), respectively 
(c.f., Bal bus&Henril l2008). The mass accretion rate is fixed 
as the Eddington one, that is m = 1 in both figures. The nor- 
malization parameter ao is an arbitrary constant and is deter- 
mined for the viscosity parameter a being unity at its max- 
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FIG. 1. — Radial profiles of the magnetic Prandtl number (left top), a 
(right top), the radial velocity (left middle), the scale height (right middle), 
the surface density (left bottom), and the temperature (right bottom). «o i- s 
determined so that the maximum of a is unity. Thick solid curves denote 
the solutions for the conventional standard disk model, while dashed, dotted, 
and dash-dotted ones do those for 6 = 0.25,0.5,1, respectively. The other 
parameters are m= 10, rh= 1, 

imum. Solid curves denote the solutions for S = (conven- 
tional standard a-model), while the dashed, dotted, and dash- 
dotted curves are for the cases S = 0.25 , 0.5 , 1 .0 respectively. 

We find that the magnetic Prandtl number is smaller than 
unity in the outer part of the disk, while it increases with 
decreasing the radius. Around the middle part of the disk, 
there is a critical radius r c where the magnetic Prandtl num- 
ber exceeds unity. It exists around r c ~ 100r s for both mod- 
el s. This is qualitatively consistent with the result obtained 
in lBalbus & Henril (1200 8). The kinematic viscosity increases 
with the disk temperature and the magnetic diffusivity de- 
creases with it. Since the disk temperature is a decreas- 
ing function of the disk radius, the magnetic Prandtl num- 
ber varies sharply with the disk temperature and becomes 
lower and lower with the increase in the disk radius [see equa- 
tion (01] except the inner disk. It would be important to note 
that the magnetic Prandtl number is insensitive to the index 5 
although the disk temperature largely changes depending on 
it in the outer part of the disk. This is because the surface 
density varies more rapidly with the disk radius compared to 
the temperature not for violating the energy conservation law 
[see equation ([8)]. 

The efficiency of the angular momentum transport and the 
viscous heating, which can be represented by the viscosity pa- 
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FIG. 2. — Radial profiles of the magnetic Prandtl number (left top), a 
(right top), the radial velocity (left middle), the scale height (right middle), 
the surface density (left bottom), and the temperature (right bottom). c*o is 
determined so that the maximum of a is unity. Thick solid curves denote 
the solutions for the conventional standard disk model (S = 0), while dashed, 
dotted, and dash-dotted ones do those for <5 = 0.25,0.5, 1, respectively. The 



other parameters are m = 10 



1. 



rameter a, must be lower in the outer part of the disk in our 
model. In addition, the parameter a strongly depends on the 
index 5. We stress here that, for the model with the larger in- 
dex S, the viscosity parameter a changes more steeply with 
radius. This tendency can be understood from equation (fl3 
by considering that magnetic Prandtl number is insensitive to 
the index S. The lower transport efficiency of the angular mo- 
mentum results in the higher density and the lower accretion 
velocity in the outer part of the disk to ensure the constant 
mass accretion rate [see equation @]. The drastic structural 
change in the middle and outer regions of the disk is a natural 
consequence of the Pm-dependence of the viscosity parameter 
a. 

We mention the physical behaviors of the inner disk. In our 
model, there exists a maximum of the parameter a at around 
the innermost radius R m . The disk structure is drastically 
changed across this point. Since we ado pted the zero stress 
bound ary condition at the inner radius dNovikov & Thornel 
fl973l) . the gaseous material cannot rotate with the Keplerian 
velocity inside R[ n . The gas would thus accrete to the central 
object while being less affected by the turbulent viscous heat- 
ing in such a region. This results in the low disk temperature 
and then the low transport efficiency of the angular momen- 



tum (plunging region). This should be the reason why there is 
an inflection point in the disk structure in our model. 

3. STABILITY OF DISK MODEL 

The inclusion of the Pm-dependence in the viscosity param- 
eter a makes drastic changes in the disk structure. We investi- 
gate the linear stability of the disk model with "non-standard" 
a-prescription to the gravitational, thermal, and secular insta- 
bilities. 

3.1. Gravitational Instability 

As is shown in § 2, the surface density in our disk model is 
larger for a larger S, especially in the outer region, than that in 
the standard a-disk model. This is because the transport effi- 
ciency of the angular momentum becomes lower in the outer 
disk when we consider the /^-dependence of the viscosity 
parameter. The larger amount of the gaseous material must be 
distributed there for the larger value of the index 6. It would 
be thus affected more by the self-gravity in our disk model 
than the standard a-model. 

We restrict our attention to the disk whose potential is dom- 
inated by the central object and whose rotation curve is there- 
fore Keplerian. The gravitational instability to axisymmetric 
perturbations then sets in when the sound speed c s , the rota- 
tion frequency 57, and the surface density £ satisfy 



Q- 



< 1 



(13) 



dToomrell96 4; Gold reich & Lvnden-Belll 1965b . The instabil- 
ity condition ( fT3l can be rewritten, for a disk with scale height 
H ~ c s /Vt around a central object of mass M, 



M disk > —M , 
r 



(14) 



where Md; s k = 7rr 2 E (see iJohnson & Ga mmie 2003, for de- 
tails). The disk becomes unstable to the gravitational instabil- 
ity when the self-gravity of the gaseous material overcomes 
the gravitational force provided by the central object acting 
upon it. 

Figures [3] and @] show the radial profiles of Toomre's 
Q-value for the models m = 10 and m = 10 8 . Solid, 
dashed, dotted, and dash-dotted curves denote the cases S = 
0,0.25,0.5, 1.0, respectively. The mass accretion rate is fixed 
as m = 1 in both figures. It is found from these figures that 
Toomre's Q-value is smaller for the model with the larger in- 
dex (5. This is because the larger amount of gaseous material 
with the slower accretion velocity should be in the outer re- 
gions for the larger index 5. Hence, the disk becomes more 
unstable to the gravitational instability for the models with the 
larger index S. 

When m = 10, Toomre's g- va hi e becomes smaller than 
unity in the outer part where r > 10 8 r s . The X-ray bina- 
ries, which are the binary system of the stellar mass black 
hole and the companion stellar object, would be gravitation- 
ally stable because the typical size of these systems would be 
much smaller than 10 s r s . In contrast, for the model m = 10 8 
(AGNs), Toomre's Q-value is smaller than unity in the region 
r > 10 2 -10 3 r s . This indicates that the geometrically thin disk 
characterizing AGN systems is gravitationally unstable, not 
only in our model, but also in the conventional a model. The 
Pm-dependence of the viscosity parameter, which is the as- 
sumption based on the recent MRI studies, enhances the self- 
gravitational instability of the accretion disk because of the 
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FIG. 3. — Radial profile of Toomre's Q-value for m = 10. Solid, dashed, 
dotted, and dash-dotted curves represent for 5 = 0, 0.25, 0.5, 1 .0, respectively. 

lower transport efficiency of the angular momentum in the re- 
gion characterized by the lower magnetic Prandtl number. 

The self-gravity should play a significant role in enhanc- 
ing the tu rbulent viscosity, and triggering the disk frag- 
mentation (Bodenh eimer et al.1 fl980t IShore & White! Il982t 
iLin & Pringldll987tlShTosman & Begelmanlll987l) . Although 
we must take account of these self-gravity effects consis- 
tently in our governing equations to construc t the disk model 
corre c tly capturing the AGNs disk system ( Kozlows ki et al.1 
Il979t [Mineshig e~& Umemuralll996|) . it is beyond the scope 
of this paper. The main purpose of this paper is to extend 
the standard a-disk to the model with Pm- dependence of the 
viscosity parameter as the first step. 

3.2. Thermal and Secular Instabilities 

Equations ([TJi and indicate that the lower disk tem- 
perature provides the lower transport efficiency of the angu- 
lar momentum, and therefore the lower viscous heating rate. 
The radial variations of the viscous parameter associated with 
radially declining viscous heating would impact on the re- 
sponses of_fhe_dj^k_jrrodd_to the thermal and secular insta- 
bilities dBalbus & Henrill2008l) . 

Since the thermal timescale is much longer than the dynam- 
ical one, we can assume that the dynamical equilibrium is re- 
tained when we consider the thermal instability. In addition, 
we can assume that the surface density of the disk does not 
change in the thermal evolution time because the centrifugal 
force balances with gravitational force in the radial direction. 
Thus the gaseous motion caused by the temperature pertur- 
bation is restricted almost in the vertical direction. These as- 
sumptions can simplify the treatment of the thermal instabil- 
ity. 

We consider the thermal equilibrium as the unperturbed 
state in the condition with the fixed surface density. When the 
disk temperature is slightly perturbed over that of the equilib- 
rium state, the criterion for the thermal instability is given 



dT 



>0. 



£=const 



where 



rad " 



32oT 4 
3r 



(15) 

(16) 
(17) 
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FIG. 4. — Radial profile of Toomre's Q-vaiue for m = 10 s . Solid, dashed, 
dotted, and dash-dotted curves represent for j = 0. . 25.0.5. 1 .0. resp ectively. 

(Shak ura & Sunvaevlll973t lPringld[l976l lPkanlll97 j). Here 
<2vi s is the viscous heating rate (Pringlel fl976l) . and Q~ ad is 
the radiative cooling rate. The thermal equilibrium satisfies 
the condition Q* is = Q~ d . Note that the viscous heating and 
radiative cooling rates are both the increasing function of the 
disk temperature. The disk becomes thermally unstable when 
the increasing rate of the viscous heating overcomes that of 
the radiative cooling. Once the disk temperature is slightly 
increased from the equilibrium state, it is further enhanced by 
the positive feedback from the visc ous heating which is more 
effective than the radiative cooling dKato et al.lll998h . 

The secular instability is a phenomenon resulting from the 
spatial modulation of the accretion rate. The typical timescale 
in which the instability grows is thus the viscous timescale, 
which is much longer than the thermal and dynamical ones in 
the geometrically thin disk. Thus we can assume that the state 
is in thermal and dynamical equilibria, so that the thermal and 
oscillatory modes can be filtered out. Then the criterion for 
the secular instability is 



dT r4 



as 



<0. 



(18) 



e:, s =e ra 



dLightman & Eardlevl [19741: iKato et alj fl998h where T r<t> = 
-2apH is the r</>-component of the turbulent stress tensor. 
This means that the instability sets in when the turbulent stress 
decreases with increasing the surface density. Once the sur- 
face density is perturbed positively from the dynamical equi- 
librium state, the turbulent stress and the efficiency of the an- 
gular momentum transport decreases. The perturbation of the 
surface density is thus further enhanced by the positive feed- 
back (~ negative diffusion). 

For the simple treatment of two types of instabilities above, 
we introduce dimensionless parameters f3 and 7. The ratio of 
the gas pressure p„. ds to the total pressure p, which is the sum 
of the gas pressu re and the ra diation pressure p ra d, can provide 
the parameter f3 (Piran 1978) 



P 



(19) 



The r-(j> component of the stress tensor T rc f, has a general form 



-laplzsP 1 ' 



(20) 



From these definition, the parameters j3 and 7 range < 
ft, J < 1- The solution of 7 = corresponds to that of the 
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FIG. 5. — Thermal and secular unstable domain on the parameter space of 
0—"f — 5. We assume electron-electron scattering for the main opacity. The 
shaded domain shows the stable region. 

Shakura-Sunyaev disk, and 7 = 1 to that where T r( p is propor- 
tional to the gas pressure instead of the total pressure. 

By adding the small-amplitude perturbations to the sys- 
tems described by set of equations d4li-([T2li including the in- 
dependence of the viscosity parameter a with the assumptions 
discussed above, we obtain the criteria for the thermal and 
secular instabilities with the parameters 6, (3, and 7, (see, ap- 
pendix for derivation). 

When we assume the electron scattering mainly contributes 
to the opacity, we can combine the criteri a for the thermal a nd 
secular instabilities to a single condition dKato et al|l9 98) as 



4-10/3-77(l-/3) + <5(8 + £)>0. 



(21) 



Note that the a-prescription depending on Pm is adopted here 
instead of the conventional approach. Neglecting the last term 
resulting from the Pra-dependence of the viscous parameter, 
equation ( f2TT > reduces to the criterion of thermal and secular 
instabilities for the standard a-model. Since the last term is 
always positive for the case 5 > 0, our disk model with 'non- 
standard' a-parameter should be more destabilized when the 
transport efficiency of the angular momentum by the MRI- 
driven turbulence depends more strongly on the magnetic 
Prandtl number. 

Figure [5] shows the stability criterion for the thermal and 
secular instabilities in the parameter space of j3—j—S. The 
shaded domain indicates the stable region. The disk is unsta- 
ble for the model with the smaller (3 and 7. This indicates 
that the instabilities are generally facilitated when the contri- 
bution of the radiation pressure to the total pressure becomes 
larger. It is characteristic of our disk model that the unstable 
domain extends for a larger 6, indicating that the disk is more 
unstable when the magnetic Prandtl number strongly affects 
the angular momentum transport in MPJ-driven turbulence. 

Considering the Shakura-Sunyaev type stress tenor (7 = 0), 
the condition for the instability is rewritten by 

The disk becomes unstable when f3 < 0.4 in the condition S = 
(standard a-model). The unstable parameter space extends 
when S > 0. This can be understood intuitively from equa- 
tions (Q~|) and (01). Our model gives the relation a cx p~ s T AS 
and thus the viscous heating rate increases with the disk tem- 
perature when 6 > 0. Since the larger index 6 provides more 
efficient viscous heating, it should play a significant role in en- 
hancing the thermal instability. On the other hand, the trans- 
port efficiency of the angular momentum decreases with the 
density when S > 0. The positive feedback effect becomes 
larger for the model with the larger 6. Hence the increase of 
the index 6 also enhances the secular instability. 




FIG. 6. — Thermal and secular unstable domain on the parameter space of 
— 7 — 8. We assume free-free opacity. The shaded domain shows the stable 
region. 

We would like to stress here that the disk becomes unstable 
even if the radiation pressure is negligibly small (fi = 1) when 
the condition S > 2/3 is satisfied. Recent local shearing box 
simulations for the MRI imply that the index 5 ranges from 
0.25 to unity dLesur & LongarettH 120071 : ISimon & Hawlevl 
2009). This suggests that the accretion disk sustained by the 
MRI-driven turbulence can become thermally and secularly 
unstable even when the gas pressure dominates the radiation 
pressure (/? gas 3> /? ra d). This is an important property that does 
not appear in the standard a-model. 

When the free-free emission is the dominant source of the 
opacity, the criteria for the thermal and secular instabilities 
are combined again to a single condition (see appendix for 
derivation) 



7 + 21/3+ 14-7(1 -/?)- 2(5(8 + /3) < , 



(23) 



Note that all terms except the last one take always positive. 
This suggests that the disk is always stable to the thermal and 
secular instabilities in the condition S = when the free-free 
opacity becomes dominant. In contrast, the last term is always 
negative when 5 > 0. The outer region of our model then tends 
to be unstable to these instabilities even when the opacity of 
the disk is mainly contributed from the free-free absorption. 

Figure [6] demonstrates the stability criterion for the thermal 
and secular instabilities in the parameter space of j3—y-8. The 
shaded region indicates the stable domain. It is found that 
the parameter space for the disk being unstable to these in- 
stabilities extends with increasing the index 6. It is important 
that the disk becomes unstable to these instabilities in the case 
5 > although it is thermally and secularly stable in the whole 
parameter domain when 6 = 0. 

The Shakura-Sunyaev type stress tensor (7 = 0) reduces the 
condition ( f23l to 

The conventional a-disk (S = 0) is always stable because f3 
should be positive. In the case 6 = 0, the radiative cooling 
rate Q~ ad more steeply responds to the disk temperature than 
the viscous heating rate g+ is . The thermal instability is thus 
suppressed by the negative feedback resulting from the strong 
radiative cooling. When the viscosity parameter a depends 
on the magnetic Prandtl number, the viscous heating rate also 
becomes the function of the disk radius. Then the positive 
feedback effect is amplified and the disk becomes unstable 
because the viscous heating rate increases with the disk tem- 
perature more steeply than the radiative cooling rate. Since 
the denominator of equation d24T i has a positive value accord- 
ing to the recent local shearing box simulations, the condition 
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d24b reduces to f3 < 16(5-7. Then the disk can be thermally 
and secularly unstable in the case S > 7/16 when the radia- 
tion pressure dominates the gas pressure (/? = 0) with Shakura- 
Sunyaev type stress tensor (7 = 0). 

When the gas pressure dominates the radiation pressure 
(f3 = 1), the criterion for the thermal and secular instabilities 
is given as 8 > 14/9 from equation d23l l. Thus the outer part 
of the accretion disks can be unstable to these instabilities 
when the efficiency of the angular momentum transport and 
the viscous heating rate are strongly controlled by the mag- 
netic Prandtl number. 

4. SUMMARY AND DISCUSSION 

We construct a disk model with the non-standard a- 
prescription depending on the magnetic Prandtl number a oc 
Pm s according to the recent numerical studies for the MRI. 
The magnetic Prandtl number Pm evaluated from the Spitzer's 
value is a decreasing function with disk radiu s and it becomes 
small er than unity for the region r > I00r s (Bal bus & Henril 
l20Q8h . 

As a result of /^-dependence of the viscosity parameter 
a, the transport efficiency of the angular momentum becomes 
lower in the outer part of the disk in the model with larger 
value of the index 8. Since the accretion velocity decreases 
with 8 due to the inefficient angular momentum transport, 
there should be a large amount of gaseous materials in the 
outer part of the disk. The self-gravity enhanced as a re- 
sult of the Pm-dependence of a-parameter would make the 
disk more unstable to the self-gravitational instability espe- 
cially for the AGN disks in the region r > 10 2 — lO 3 ^ (see 
Paczynski|[l978t iKozlowski et al.lll979t iFukue & Sakamoto! 
1992t Mineshig e & Umemural 119961" for the case with stan- 
dard a-prescription). 

Since our model gives the relation a oc p~ s T 4S , the viscous 
heating rate increases with the disk temperature when 8 > 0. 
This suggests that the positive temperature perturbation en- 
hances the viscous heating rate and thus the positive feedback 
to the thermal instability. On the other hand, the transport ef- 
ficiency of the angular momentum decreases with increasing 
the gas density when 8 > 0. This works in the same way as 
the negative diffusion and also enhances the positive feedback 
to the secular instability. 

It is the most remarkable feature of our disk model that the 
thermal and secular instabilities can grow in its middle part 
even if the radiation pressure is negligibly small in the condi- 
tion 8 > 2/3. Also the outer part of the disk can be unstable 
to thermal and secular instabilities when 5 > 7/16 (radiation 
pressure dominant with Shakura-Sunyaev type stress tensor, 
7 = 0) or 8 > 14/9 (gas pressure dominant). If the MRI plays 
an important role in operating the MHD-turbulence and the 
viscosity parameter a depends on the magnetic Prandtl num- 
ber, the steady mass accretion would not be maintained due 
to the growth of these instabilities in the geometrically thin 
disk. The rapid transition from the steady accretion phase to 
the non-steady state should be considered to enable us to un- 
derstand more deeply the origin of the light curve observed 
from the realistic astrophysical disk system. 

In our model, we assumed that the turbulent stress tensor 
is just proportional to the pressure according to the standard 
a-disk model (see, equation l20l . On the other hand, re- 
cent local shearing box simulations for the MRI imply that 
the pressure dependence of the turbulent viscosity might be 
weak compared to that assumed in the standard a-disk model, 
that is T n< p oc p e where the index e takes a value in the range 



1/6 < e< 1/4 (ISano et al.ll2004t iPessah et al.ll2007l) . 

When the pressure dependence of the turbulent viscosity is 
weaker than that assumed in the standard a-disk, the viscous 
heating rate is expected to become insensitive to the temper- 
ature perturbation. The disk thus becomes more stable to the 
thermal instability. Additionally, the secular instability should 
also be suppressed because the turbulent stress would then re- 
spond weakly to the density perturbation. 

Now let us derive the stability condition of disks when the 
turbulent stress varies weakly depending on the pressure as is 
implied from the recent numerical studies. Using the prescrip- 
tion T r< j, = -2af>gas/? e H [note that this equation reduces equa- 
tion (|20| > when e = 1 -7] and assuming the electron scatter- 
ing as the dominant opacity source, the instability conditions 
are given by ( IA7t for the thermal instability, and (IA121 > for 
the secular instability (see Appendix for the details of deriva- 
tion). Conditions iA7l and JA121 > can commonly yield, by 
taking (3=1 (gas pressure dominant), 

7 — 7 — e 



S> 



(25) 



Here we assume 8 > 0. This equation indicates that the weak 
dependence of the viscosity parameter on the pressure (i.e., 
I7I, |e| <C 1) has a stabilization effect on the thermal and sec- 
ular instabilities. 

For the model in which the free-free absorption is the main 
opacity source, the stability conditions for the thermal and 
secular instabilities can be obtained from JA9t and ( 1A131 >. re- 
spectively. By taking (3 = 1, these conditions are combined to 
a single condition 

S> 1 -^. (26) 

Here we assume 8 > 0. This also shows that the instabilities 
can be suppressed when the pressure dependence of the vis- 
cosity parameter is weaker than that assumed in the standard 
disk models. 

It should be important to stress that even when the turbulent 
stress does not directly depend on the pressure (i.e., 7 = e = 0), 
the disk can be unstable since the efficiency of the angular 
momentum transport a depends on the temperature and the 
density. In such case, the instabilities grow when 8 > 7/9 for 
the middle part of the disks and when 8 > 5/3 for the outer 
part of the disks. Since 8 is predicted to be less than unity from 
the recent local simulations, the middle part of the disks can 
be unstable to these instabilities when the saturation levels of 
MRI-driven turbulence is strongly controlled by the magnetic 
Prandtl number. 

Finally we discuss the validity of the assumption adopted 
in our disk model. Numerical studies for the MRI imply that 
transport properties of the MRI-driven turbulence depend on 
the magnetic Prandtl number as is given in equation (HJ. 

This relation is verified only in the narr ow parameter range 
10" 2 5, Pm < 10 and 10 < Am < 10 4 dLesur & Longarettil 
I2007t lLongaretti & Lesuri |2010h . Here R M = Hcjr] is the 
magnetic Reynolds number. The range of Pm examined in 
numerical studies corresponds to the region 10r, < r < 10 3 r, 
in our model. The fm-dependence of the viscosity parameter 
a is not well studied due to the computational issue in the pa- 
rameter range characterizing the outer part of the disk model. 
In addition, Reynolds and magnetic Reynolds numbers are 
taken to be involuntarily larger values than the realistic val- 
ues for the astrophysical ionized gas due to the limitation of 
the computational resources. If we suppose the realistic val- 
ues of viscosity and resistivity, their related dissipation scales 
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become much smaller than those resolved in existing state- 
of-the-art numerical simulations. In order to facilitate the un- 
derstanding of astrophysical accretion processes, we need to 
examine whether such small scale dissipation processes actu- 
ally affect the dynamics of the accretion disks via the angular 
momentum transport and then the instabilities of the disk. 

Additionally, it has not been confirmed whether the vis- 
cosity parameter is impacted by the magnetic Prandtl num- 
ber when the microscopic diffusivities are the functions of 
the physical parameters such as the density and the temper- 
ature because the uniform diffusivities are assumed in the 
existing local shearing box simulations (|Les ur & Lo ngarettH 
120071: iFromang et alJ2007tlSimon & Hawleyl2009l) . We need 
further systematic studies on the non-linear properties of the 
MRI in the vertically stratified disks with variable diffusive 
parameters and the realistic radiative cooling in order to ver- 
ify the Pm-dependence of the turbulent viscosity parameter 
a. 



Our model does not include just for simplicity the effect 
of magnetic fields explicitly (but takes account of the ef- 
fect of the MHD turbulence through the non-standard a- 
prescription). The impacts of the magnetic field can modify 
the disk structure. Since the magnetic pressure contributes 
to the total pressure appeared in equation d20T ) and weakens 
the temperature dependence of the viscous heating, it might 
suppress the growth of gravitational, thermal and secular in- 
stabilities when the magnetic pressure dominates the gas and 
radiation pressures. This is beyond the scope of this paper but 
would be studied as our future work using numerical proce- 
dures. 



We are grateful to an anonymous referee for improving our 
manuscript. We also thank Hiroshi Oda, Ken Ohsuga and Ry- 
oji Matsumoto for fruitful discussions. 



APPENDIX 
CRITERIA FOR THE INSTABILITIES 



In this appendix, we briefly summarize the derivation of the stability conditions shown by equations ( l2Tb . d23l l. ( f25l l. and ( f26*b . 
Here we use a following general prescription for the stress tensor, instead of the form introduced by equation (F20l i. 



T r4> = -2ap ] p £ H , 



(Al) 



where e is a constant. When e = 1 -7, this equation reduces to equation d20b . We hereafter describe the unperturbed and perturbed 
quantities by subscripts '0' and '1', respectively. Note that all the physical unperturbed quantities are assumed to be much larger 
than their associated perturbed quantities in the following. 

We consider the stability criterion of the system to the thermal instability. As is discussed in § 13.21 we can assume that the 
dynamical equilibrium is maintained during the growth of the thermal instability since the dynamical time scale is much shorter 
than the thermal one. Also we can assume that the material motion resulting from a temperature variation occurs only in vertical 
direction since the force balances with the gravitational force in the radial direction. Then the surface density S is approximately 
unchanged during the growth of the thermal instability; 



S = const, 



while the hydrostatic balance is maintained (Kat o et aTlll998h . 
The linearized form of equation (O and the total pressure p = p s - ds + /? ra( j, where /? ra d is the radiation pressure, are 



El 

PO 



1-/3 Si 4-3/32] 
'T+~0% + 1+0 ¥ : 

20 Si 4-3/3 T { 



1+/3S 1 + T 



(A2) 



(A3) 



(A4) 



Note that equations (fl2l) and ( fl9b are adopted for the derivation. 
Then the perturbative equation of the viscous heating (fl6] i can be obtained from £[), (IAU . and ( IA3b - (IA4b as 



Q 



vis.l 



(4 -3/3)(l + e)- 7(3 -40) + 5(&+0)Ti -1 + 0-25+2^+206^! 



Q + vi S> o l+£ 7b 1 + So' 

When the electron scattering is the dominant source of the opacity, the radiative cooling (TTTb provides a linearized equation 

Si 
S ' 



(A5) 



6 ra d.l _ qT\_ 
2rad,0 T 



(A6) 



By substituting equations jA5i and jA6i into equation ( fT3T > and using equation ( 1A2I >, we obtain the instability condition 

4 - 1 0/3 - 7 7 ( 1 - 0) + 5(8 + 0) - (4 - 3/3)( 1 - 7- e) 



1+/3 



>0. 



(A7) 



When we take e = 1 - 7, this equation reduces to equation (l2TT i since > 0. 

When the free-free absorption mai nly d etermines the opacity of the system, we can give a linearized equation for the radiative 
cooling, by using equations (flTT i and dA3b . 

67ad,i _ 23 + 9/3 ri 3 + /3S, 
£mdo 2(1+0)7;, 1+/3S - 
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By substituting equations ( IA5t and (IA8b into equation (TT5T > and applying equation (IA2b . we obtain the instability criterion for the 
thermal instability, 

7 + 21/?+14 7 (l-/?)-2 ( 5(8 + /3) + 2(4-3/3)(l-7-e) , A 

2(1 + 13) < °- (A9) 

Note that this can reduce to the condition ( 1231 when e = 1 - 7 since /3 > 0. 

Finally we consider the stability condition of the system against the secular instability. When the secular instability grows, the 
negative diffusion amplifies the surface density so that the assumption (lA2b cannot be adopted here. As mentioned in § 13.21 the 
viscous time scale, which is equivalent with the typical growth time of the secular instability, is much longer than the thermal 
time scale. The thermal equilibrium is thus maintained during the development of the secular instabilities, i.e., Q+ ls = Q~. dd - This 
can provide a linearized equation 

(-4 + 7+4,5) J + (l+ 7 -«5)|i + (l- 7+( 5)^ + e^ =0, (A10) 

for the case with the electron scattering as the main opacity source, and 

(-15 + 27+85) J + (4 + 27-2<5)|i+2(-7+<5)^ + 2e^ =0, (All) 
^0 ^0 no po 

for the model with the free-free absorption as the dominant opacity source. 

By combining equations (1A31 >. (lA4b and (lAlOt or (IA1 j} , we can express the quantities H\/Hq, Ti/Tq and pi/po by Si/So 
for each model with different opacity source. Then by substituting these obtained equations into equation ( fT8l , we obtain the 
instability criterion for the secular instability. 

When the electron scattering is the main opacity source, the instability condition is given as 

5 7 +4e+/3(l+4 7 +5e+*) >Q 



4-10 / S-7 7 (l-/3) + <5(8 + £)-(4-3/3)(l- 7 -e) 

This condition is identical with equation JA7b when 5 > since the numerator in equation (IA12b has a positive value. Especially 
when e = 1 -7, this equation reduces to equation (fjj}. 
The disks in which the free-free absorption becomes the dominant opacity source is unstable to the secular instability when 

l+28 7 +24e + 26+£(3 + 8 7 +12e + 2<5) 
7 + 21/3+ 14 7 (l-/?)-2<5(8 + /3) + 2(4-3/3)(l- 7 -e) < 

Since the numerator of equation jAl 3I > is positive when 5 > 0, the condition for the secular instability in the outer region is 
identical with that for the thermal instability given in equation ( |A9t . When e = 1 -7, this equation reduces to equation ( f23l . 
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